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The Lanczos-Lovelock models of gravity constitute the most general theories of gravity in 
D— dimensions which satisfy (a) the principle of of equivalence, (b) the principle of general co- 
variance, and (c) have field equations involving derivatives of the metric tensor only up to second 
order. The mth order Lanczos-Lovelock Lagrangian is a polynomial of degree m in the curvature 
tensor. The field equations resulting from it become trivial in the critical dimension D = 2m and 
the action itself can be written as the integral of an exterior derivative of an expression involving 
the vierbeins, in the differential form language. While these results are well known, there is some 
controversy in the literature as to whether the Lanczos-Lovelock Lagrangian itself can be expressed 
as a total divergence of quantities built only from the metric and its derivatives (without using the 
vierbeins) in D = 2m. We settle this issue by showing that this is indeed possible and provide an al- 
, gorithm for its construction. In particular, we demonstrate that, in two dimensions, R^—g = djR? 

for a doublet of functions Pj = (R° , R 1 ) which depends only on the metric and its first derivatives. 
We explicitly construct families of such R?s in two dimensions. We also address related questions 
regarding the Gauss-Bonnet Lagrangian in D = 4. Finally, we demonstrate the relation between 
the Chern-Simons form and the mth order Lanczos-Lovelock Lagrangian. 
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I. INTRODUCTION AND MOTIVATION 

Consider a theory of gravity in D dimensions based on the action 



A = I V=9d D x[L(R a bcd ,g ab ) + L matt (9 ab ,M}, (1) 

where L mat t(g ab , 4>a) is a suitable matter Lagrangian involving some matter degrees of freedom &a ■ Varying g ab in 
this action, with suitable boundary conditions, leads to the field equations (see e.g. sec 3.5 of ref. [l| or, for a textbook 



description, chapter 15 of ref.[2|): 



Gab — P a Cde Rbcde — ^^dab — 2\/ € \7 d P acdb = 1Z a b — ^Lg a b — 2V C V d P ac db — -^T a b, (2) 



where 



pabcd = _ 7 ^_ (3) 
OKabcd 

We use a signature (—,+,+,...) and units with Gn = c = 1; Latin indices run over (0, 1, 2, D— 1). The notation in 
terms of calligraphic fonts is suggested by the fact that in Einstein's theory, with L = R/IQ-k, we have lZ a b = i? a &/167r 
and Gab = Gab/lQn with standard normalization. The field equation in Eq. ^ will not contain terms involving 
derivatives of the metric of degree greater than two if V a P abcd = 0. Since P abcd has the symmetries of the curvature 
tensor, it follows that it must be divergence-free in all indices. One can show that the most general tensor satisfying 
such conditions, built from the metric and curvature tensor, must be a polynomial with the indices of curvature 
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tensors contracted out in a specific manner using the determinant tensor. In such a case, the gravitational part of the 
Lagrangian will take the form 



K 



L = J2 ; l m « ^X::Cr<a 2 • • • > ( 4 ) 

m— 1 

where 8%^'"£* m is the determinant tensor (for example, ^ = ^(S^d, ~ ^d^c)): tne c «i are arbitrary constants and 
Lf m ) is called the m-th order Lanczos-Lovelock Lagrangian [3]. The m = 1 term is proportional to S^R^ oc R and 
leads to Einstein's theory. The m = 2 term gives rise to what is known as Gauss-Bonnet theory with 

L GB = L 2 oc (R 2 - 4R ab R ab + R abcd R abcd ). (5) 

In the language of differential forms and tetrads, denoting our cotangent basis w a , spin connection tt AB , and curvature 
two- form Q AB , and defining e 0l ... 0n = *(w 01 A- • • Aw 0n ), the mth Lanczos-Lovelock Lagrangian form C m can be defined 
as: 

C m = & AlBl A • • • A 9 A '" S '" A eA lBl ...A m B m . (6) 

(See, for example, [|J; we will try to consistently use calligraphic font for the differential form corresponding to the 
Lagrangian density and normal font for Lagrangian which is scalar; in most case, the context will make clear what 
is the expression we are talking about). Lanczos-Lovelock Lagrangians are the most natural generalizations of the 
Einstcin-Hilbert Lagrangian to higher dimensions, and both share several remarkable properties. It has been shown 
Q that these Lagrangians embody a version of the holographic principle, at the level of the action, and the resulting 
models share many of the thermodynamic features of the standard General Relativity @, 0] • 

Because of the determinant tensor in Eq. Q , it is obvious that in any given dimension D we can only have K terms 
where D > 2K . We will call D = 2m the critical dimension for the mth order Lanczos-Lovelock Lagrangian. It is 
possible to show by straightforward combinatorics (see e.g. that Gab = identically in the critical dimension; this 
generalizes the well known result that G a b = identically in D — 2, which is the critical dimension for Einstein- Hilbert 
action corresponding to the m = 1 Lanczos-Lovelock Lagrangian. It follows that the field equations are vacuous for 
the Lanczos-Lovelock theory at the critical dimension. 

This raises the question: What is the structure of Lanczos-Lovelock action itself in the critical dimension? The 
vacuous nature of the field equations tells us that the variation of the action will be a pure surface term in the critical 
dimension and can be ignored (or compensated) under suitable boundary conditions. But it is not obvious whether 
the action itself is a pure surface term. This is the central question we address in this paper. We begin with a series 
of comments to convince the reader why this issue is interesting and nontrivial. 

Let us focus our attention on D — 2, which is the critical dimension for m — 1 Lanczos-Lovelock theory, the 
Lagrangian of which is simply R and the action is just the standard Einstein-Hilbert action. It is trivial to show, 
using the language of differential forms, that the integrand of the Einstein-Hilbert action R^/—g in D = 2 can be 
written as an exterior derivative of a quantity which explicitly involves the tetrads (see Eq. ([8]) below). But what 
we would like to know is whether the integrand of the Einstein-Hilbert action R^J —g can be expressed as djB? for 
a doublet of functions (R r,R ) built only from the metric and its derivatives without using tetrads in a general 
coordinate system. That is, if we take the line element to be 

ds 2 = A(t, x)dt 2 + 2C(t, x)dxdt + B(t, x)dx 2 , (7) 

(in which we have not imposed any gauge or coordinate conditions) we want an expression for (i? ,^ 1 ) in terms of 
A, B, C and their derivatives. This is again trivial to accomplish if we use the conformally flat gauge with C = 0, A = 
B] then one gets R^f—g — djd J \nA which has the required form but this will not survive when we transform to an 
arbitrary coordinate system. 

Given the fact that the answer in the language of differential forms is already known, one would have thought 
that the expression for R? would already exist in the literature. We found, to our surprise, that this is not the case! 
Further, a paper by Deser and Jackiw [8| seems to suggest that this is impossible, which had led Kiriushcheva and 
Kuzmin [10| to recently conclude that there were fundamental differences between the vierbein formalism and the 
standard formalism in General Relativity (most notably, the vierbein formalism is known to have an extra internal 
gauge symmetry) and, because of those differences, the Einstein-Hilbert Lagrangian being a total derivative in the 
vierbein formalism does not imply that it is a total derivative in the standard formalism. It is therefore clear that 
this issue deserves clarification: on one hand, several papers claim (at least implicitly) that the action itself is a total 
derivative for all Lanczos-Lovelock models in critical dimensions but without concrete proof; on the other hand, there 
are arguments, summarized in the articles by Kiriushcheva and Kuzmin [9l.[l0|. that these claims are ill-founded. This 
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confusion is in large part due to the lack of an explicit expression for R 3 , as a function of the metric and its derivatives, 
satisfying L m ^J—g = djR 3 . One of our conclusions is that R 3 can indeed be expressed in terms of the metric and its 
derivatives. We will provide an explicit expression for R 3 thereby putting an end to the debate in the case of D = 2. 
This R 3 is highly non- unique: we can always add to it any set of functions f J (x) which satisfy djf 3 = 0. We will see 
that there is one natural R 1 associated with each unit normalized vector field in the D = 2 spacetime and, of course, 
any two of them differ by an f 3 with djf 3 — 0. 

Let us next turn to the case of D = 4. The Lanczos-Lovelock Lagrangian which becomes critical in D = 4 is 
the m = 2 Gauss-Bonnet term, given by Eq. ([S]). The question again is whether one can write Lgb \/ — = djR 3 . 
Such a claim, with an expression for R 3 , has been given in [ll| and has been cited in a few later papers as well 
based on an expression originally given in [l2l |. Unfortunately, this expression is incorrect and arose from an invalid 
identification of frame indices with spacetime indices. We will clarify this issue as well and will provide a proof that 
every critical-dimension Lanczos-Lovelock Lagrangian can be written as a total derivative of functions of the metric 
and its derivatives, and that these functions can be thought of as the coordinate and local representation of the 
familiar Chern-Simons form. This will be done in a pedagogical manner to ensure that the issues are clarified. 

We will now proceed to establishing these claims. The plan of the paper is as follows: we begin by addressing 
the issues related to the two-dimensional Einstein-Hilbert Lagrangian in Section |TT] by explicitly finding the pair 
of functions R 3 (g,dg) which satisfy R^—g = djR 3 . We then briefly consider the four-dimensional Gauss-Bonnet 
Lagrangian in Section IIIII before moving on to the general 2TO-dimensional Lanczos-Lovelock Lagrangians in Section 
HVl Appendix lAl gives a short review of the Cartan formalism which is used throughout the paper, while Appendices 
[61 through lEl present in detail the derivations to some of our claims. 



II. EINSTEIN-HILBERT ACTION IN D = 2 



We begin by considering the two-dimensional Einstein-Hilbert Lagrangian and will obtain an explicit expression for 
R 3 , dependent only on the metric and its derivatives, such that R^J—g = djR 3 . This will provide a direct proof that 
the Lagrangian is indeed a total derivative. We will do so in two ways to provide two perspectives on how the form of 
R 3 could have been obtained: the first procedure derives R 3 through the Cartan formalism which uses the language of 
differential forms. This is quite powerful and immediately generalizes to arbitrary D. The second procedure, specific 
to D = 2, constructs one R 3 from each unit normalized vector field. Finally, we mention a procedure based on 
conformal transformation, leaving the details to Appendix [Cj These various approaches are discussed because there 
has been some controversy regarding differences between the Cartan formalism and the standard formalism; hence, it 
is useful to derive Eq. (J9j) without resorting to the Cartan formalism at all. 



A. First Method: Cartan Formalism 



If we define an orthonormal basis co a = Lofdx 1 satisfying g$j — ujfu! b rj a i, as a basis for our cotangent space, then we 
notice that the spin connection, defined by the first Cartan structural equation (see Eq. (jAll) of Appendix A) must 
be antisymmetric: Sloi = — fiio- Then, the Einstein-Hilbert two-form can be written (up to an overall sign due to the 
definition of d 2 x): 

Ry/^gd 2 x = e ab A *{uj a A u> b ) 

2 W 

= dfloid x, 

where <d ab = (dfi ab + Q a s A Q sb ) is the curvature two- form (see Appendix B for details). Therefore, finding R 3 such 
that R^J—g = djR 3 holds is straightforward: one simply needs to calculate explicitly the spin connection J7oi and 
take its exterior derivative. This is done in Appendix [B] and the result is 



R° 
R 1 



-X—goo.i + (1 - X)—gu i - 2(1 - A)g i,i + 011,0 
goo gii 



(1 - A) — g-11,0 + A — goo.o - 2A.901 

011 000 



000,1 



(9) 



which is a one-parameter family of solutions, parametrized by a constant A. This parameter should be thought of 
as characterizing the gauge freedom we have in defining an orthonormal basis for our cotangent space. The R 3 s for 
two different values of A constitute an example of non-uniqueness as mentioned in Sec. HI One can directly verify 
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that the difference p = [R 3 — R 3 (A2)] satisfies djf 3 — 0, as it should. We could have, of course, set say A = 1 
without any loss of generality using this freedom; we have retained for future convenience. (This corresponds to the 
freedom of changing Oqi to Sloi — > ^01 + / where df — 0.) It should also be emphasized that this gauge freedom is 
entirely unphysical, and one can fix it without imposing any condition on our spacetime or our metric. In this sense, 
the gauge-fixed Cartan formalism is equivalent to the standard formalism for General Relativity. 



B. Second Method: R J from normalized vectors 



In two dimensions, we know that R ab = (l/2)g ab R; this will allow us to express R as the total divergence of a vector 
field by the following procedure. Let n l be a unit normalized vector in two dimensions which we take to be timelike 
(with riiJi 1 = — 1) for definiteness; the same procedure works for spacelike vectors as well. Using the standard relation 
between n b R ab n a and the commutator of covariant derivatives, it is easy to show that 



R bd n b n d = 7i b V a V 6 n a - n 6 V fc V a ?i a = V a (n b V b n a - n a V b n b ) + M b b M a a - M\M b al 



(10) 



where we defined the tensor M ab = V a n b . The last two terms cancel out in two dimensions. To see this, note that 
[M b b M a a — M a b Al b a ] is a scalar; since n b M ab = 0, then, choosing coordinates where n a cx 6q, we find that M a0 = 
and M a = 0. This then implies 



M b b M a a - M\M b a = M\M\ - M\M\ = 0. 



Further, using R a b = (l/2).g a fc-R and niv} = — 1, we combine Eq. (fTT|) and Eq. ([TO]) to find 

R = 2V a {n a V b n b - n b V b n a ). 



(11) 



(12) 



It follows that R^ 



d a R a , where 



R a 



-g{n a V b n b -n b V b n a ), 



(13) 



for any unit normalized vector n a . Since we get one R a for each n a , there is infinite degeneracy in the expression; 
any two such R 3 will differ by an p with djf 3 = 0. 

The derivation above is almost identical to the standard Gauss-Codazzi decomposition (see e.g. (12.55) of 0) of 



R in terms of the extrinsic curvature Ki. 



-Mi, 



j-tj — Tiiin V fcTij). In two dimensions, the only non-zero component 
of is K n , which means that Tr(K 2 j = (TrK) 2 . Since we also have (TrK) 2 = (TrMf and Tr(K 2 ) = Tr(M 2 ), 
then this decomposition is identical to Eq. (fTTj) and we get R = —IVj^Kn 1 + a J ) with a J = n b V b n : > . We gave a 
direct derivation of Eq. (fT2"j) to stress that we need not introduce hypersurface orthogonality or a foliation to obtain 
the above result; the fact that M ab is not symmetric is also irrelevant. Of course, this method only works because we 
are in two dimensions. 

By working out the components of i? J with n l — (— 5oo) _1 ^ 2( ^0' one obtains 



R° = 



R 1 = 



9oi 

5oo,i + 5n,o 

5oo 



5oi , 

— 5oo,o — ^501,0 + 5oo,i 

5oo 



(14) 



which is the same as Eq. © with A = 1. Repeating the derivation with a spacelike vector ri 1 cx S\, one would obtain 
the corresponding result with A = 0. The linear combination of these solutions R 3 is also a solution, giving us Eq. (JSJ. 
It is therefore clear that R^—g can indeed be written as the total derivative of a function of the metric and its 
derivatives in any given coordinate system in an infinitely degenerate way. 

In the paper by Deser and Jackiw [8| mentioned earlier, it is stated that "....i?^ cannot be presented explicitly and 

locally in terms of the generic metric g^ and its derivatives dag^y " . One of the authors has informed us (Jackiw, 

private communication) that they were not interested in an expression for R J given in terms of components. In that 
case, the sentence quoted above does not exclude expressions such as Eq. ([§]), Eq. (fl"4"|) etc. which use the components 
explicitly. In contrast, Eq. (|13p does give an expression for R 3 explicitly as a vector density, but uses an external 
construct — the vector field n % — in addition to the metric and the connection (and hence again does not contradict 
the statement quoted above). While one cannot construct nontrivial tensorial quantities (except, of course, for the 
density \J—g) from the metric and its first derivatives, it is certainly possible to do so if additional vector fields are 
given. 
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To avoid misunderstanding we stress the following point: In any coordinate system one can choose a normal- 

ized timelike vector n l — (— goo) -1 ^ 2 ^ with components proportional to (1, 0) and use Eq. (fT3|) to obtain Eq. (fl4|) in 
that particular frame. If we now perform a coordinate transformation x % — > x 1 , Eq. (|13j) will transform in a generally 
covariant manner with R J \-\/— g)^ 1 ^ 2 transforming as a vector. However, since n l will no longer be proportional to 
(1,0), the explicit expression, Eq. (|14p . will not hold in the new coordinate system. Instead, there will exist another 
vector, N$ = (— goo) _1 ^ 2< 5o: m the new coordinate system which yields the same expression, Eq. (IT4"|) . in terms of g ab . 
(Recall that Eq. (fT2"]) holds for any normalized timelike vector field; each gives one particular B? through Eq. (p3|).) 
Similar results hold for a unit normalized spacelike vector field. In other words, in any frame, we can always choose 
a vector n l such that Eq. A13\) reduces to Eq. \1J$ . In this sense, the expression in Eq. (|14p holds in the same form in 
all frames, though in each frame it is associated with a different vector field. Therefore, expressions like Eq. ^ arc 
as good as expressions involving abstract indices (e.g. djg a b) rather than explicit components (e.g. d^goi). 

It is also possible to obtain Eq. ([9]) by demanding that B? have the correct conformal transformation properties; 
the details of this method are in Appendix [C] 



III. GAUSS-BONNET LAGRANGIAN IN D = 4 



Before moving on to the general Lanczos-Lovelock action in arbitrary dimension, we shall quickly consider the 
four-dimensional Gauss-Bonnet Lagrangian. We begin by rewriting LcB\J~^g in the differential form language: 

c GB = (dn AB + n A s a n SB )A 

{dn CD + n c T a n TD ) a e ABCD (15) 

n A B a (eg - in^) a e // 

This equation is derived (see Appendix [Dj in the orthonormal basis where the spin connection Q,ab is antisymmetric; 
it is the analogue of Eq. ([5]). Proceeding exactly as in the case of D = 2, we can, in principle, obtain an expression 
similar to Eq. ^ thereby expressing the R J in terms of the metric and its derivatives. In practice, however, this 
calculation, while conceptually simple, becomes tediously long due to the overwhelming number of terms involved. 
The reason for this is easy to understand: in the m — 1 case, we had C — dfi with only one non-zero fi a {, term: f2oi- 
This term was built from two function /j, bc a as per Eq. (|B5[) from Appendix each of which is built of up to two 
vierbein terms. We therefore end up with approximately 1x2x2 = 4 terms for B? . 

In the m — 2 case, each R? is actually built from a sum of 36 terms, each of the form VL A (dQ, + ilAft — f2 A f2) . Each 
such n is built from eight fi bc a terms, each of which is built from four terms. So at first estimate, we would end up 
with 36 x (8 x 4) 3 ps 10 6 terms. While this is a rough calculation (and overestimate), it is clear that the expressions 
will be far larger than they were in the D = 2 case. Nevertheless, we will illustrate the procedure for one particular 
metric: consider the line element 

ds 2 = g o{t)dt 2 + 2g i(t)dtdx + 2g 2(t)dtdy + 2g 03 (t)dtdz + gn(t)dx 2 + g 2 2(t)dy 2 + g 33 (t)dz 2 , (16) 

where every component depends exclusively on the time coordinate. Then, using the steps outlined in Appendix [B] 
one finds 

R° = — 322,05ll,0933,0ff00ff33ff229ll x 

-.9iig22goo + giiff^ffoo + 3ffiig22goig 2 - 2g 2 1 ff§ 2 ffg 1 .g o - 3giigg 1 g 2 2goo + .9ii.9qi.9§ 2 + ffo^ (17) 

( n n n -X- n n 2 -I- n 2 n A 5 / 2 ( n n J- n 2 \ 2 ( 900 91 1 322933 -flll gg 3 322 ~3ll 9333^2 ~g33g22 gpl \ I \ 

(-googiig22+g22.9 i+3o23ii) (-Sdoflii +5oi) { g229ooSll _ S22S ^ _ ff g aail ) {-9m) 

which satisfies 

J—g (R 2 - AR ab R ab + R abcd R abcd ) = d R°. (18) 

It would therefore be advantageous to find a way to write R^ using abstract indices, Christoffcl symbols, and other 
well-known quantities in order to make it easier to work with. Further note that — because of the non uniqueness 
of R J — it is entirely possible that the metric in Eq. (I16|) admits another R 3 , which is a lot simpler and differs from 
the one in Eq. (|17|) by a set of functions f° with djf J — 0. It is very unwieldy to start with an expression such as 
Eq. (|17p and attempt to rewrite it or simplify it. However, as we shall see in the rest of this section and in Section 
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IIV A[ there are very large classes of spacetimes which do actually allow us to write B? in rather simple forms (In fact, 
as we will see in Section IIV Al there are indeed much simpler forms of i? J available for the metric in Eq. ([T5| ) . 

Cherubini et al.[ll[ has an expression for R? for the four-dimensional Gauss-Bonnet Lagrangian which is based 
on pj)]. This expression is the only one for B? (for any Lanczos-Lovelock Lagrangian) which we could find in the 
literature: 



L 



GB\ 



-d a 



abed ranj^r 
~9V Vrs 1 mb 



2 R 



ncd 



cl rid 



(19) 



where r\ refers to the Levi-Civita symbol. We found its simplicity amazing, but we realized that — unfortunately — 
this equation is algebraically incorrect: it fails for just about any metric. By a curiously incorrect reasoning, one can 
obtain a similar expression which works for a large class of metrics [13( ■ We describe this result briefly, in view of its 
potential interest and to correct some of the errors in the formulas in the literature. 

Note that our expression in Eq. (|15p is derived in the orthonormal basis in which one cannot — in general — express 
the spin-connections in terms of Christoffel symbols by Q a b — T b l c dx c . This relation, however, is true in a coordinate 
basis (for a review of the differential form language, see e.g., chapter 11 of ref. Q or Appendix |A"|) . Interesting things 
happen if we forget this fact and substitute n a b = T bc dx c (which is certainly incorrect) for the spin connections in 
Eq. |(I5]). In that case, we find 



LGB\F r g = -d a 



-r 



' V 1 

il nd 



d a Q a - 



(20) 



The expressions in the right hand sides of Eq. (fTO)) and Eq. (l20l) differ only by a sign of a T 2 term. (This sign difference, 
it turns out, is simply due to confusion caused by differing conventions between 12 1 and [Tl|). This expression in 
Eq. (|20[) however — quite surprisingly, since we made an obviously incorrect assumption — happens to be correct for 
a large number of spacetimes, such as Schwarzschild, Kerr, Milne, Godel and even Ozvath-Schuking. This expression 
also has the advantage of being written in a concise notation: it uses abstract indices and familiar quantities such as 
Christoffel symbols and curvature tensors. This is unlike our previous formulas for BP, such as Eq. ([9]). 
But, of course, Eq. is incorrect in general: for example, a simple line element for which it fails is: 



ds 2 



g 00 (z)dt 2 + gn(t)dx 2 + gi2{y)dy 2 + 2dydz + gzs{t)dz 2 



(21) 



where we made explicit the coordinate dependence of the components. One can also use a closely related form of 
the metric to illustrate the fact that Eq. (|20]) is not generally covariant and whether it holds or not depends on the 
coordinate system used. (This is to be expected because Q a is not a vector density.) Consider the following metric, 
obtained from Eq. ((2"T]) by setting 522 = 1: 



ds 2 = goo(z)dt 2 + gn(t)dx 2 + dy 2 + 2dydz + g^(t)dz 2 



(22) 



One finds that Eq. (|20]) holds. But if we introduce a new coordinate ( by y = (l/2)£ 2 , so that the metric becomes: 



ds 2 = g 00 (z)dt 2 + gn(t)dx 2 + ( 2 d( 2 + 2(d(dz + g 33 {t)dz 2 



(23) 



then we find that Eq. (f2"0"]) does not hold! For the metrics in Eq. ((2"2"]l or Eq 
are looking for a Q- 7 with djQi — 0. For the metric in Eq. ((22]) . we find that Q- 7 
Eq. (E3]) we find 



Q° = 



-1 



911, 0000,3 



CVg, 



33 



T 3 / 2 

l9oo a 



the Lqb actually vanishes so we 
for all j. But for the metric in 



(24) 



which has dgQ ^ 0. Obviously, the expression for Q 3 is very much coordinate dependent. In any case, we do not 
expect an expression like Eq. (|2U]) to be valid because it arose through an incorrect mathematical operation. 

One of the authors of ref. 12] (in which Eq. ((20]) was first derived) has confirmed (F W.Hehl, private communication) 
that Eq. (|20l) is indeed incorrect and the error arose due to an invalid identification of frame indices with spacetime 
indices in [12]. Thus, there is no simple expression for i? J in the case of D = 4 for a general metric (though, as we 
shall show later, one can find such an expression if the metric possesses a Killing vector; see Sec. IIV Al) . 



IV. GENERAL CASE: LANCZOS-LOVELOCK LAGRANGIANS IN D = 2m 



We shall now consider the general case of Lanczos-Lovelock Lagrangians in the critical dimension D = 2m. There 
exists a simple geometrical interpretation for the Lanczos-Lovelock Lagrangians being total derivatives of expressions 
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involving tetrads. Indeed, C m is simply the Euler density of a 2m-dimensional manifold, as can be seen, for example, 
in the section on index theorems in ref. [l5j . When working in 2m dimensions, C m is a topologically invariant 
scalar density because it is a closed form and therefore locally exact [HI]. The Chern-Simons form Q m is defined to 
embody this quality: evaluating it between two connections SI and SI', one obtains C m (Q) — £ m (Sl') = dQ m (fl,tt'). 
In particular, if we restrict ourselves to a local coordinate patch, then we can introduce a flat connection SI' — and 
can therefore write C m — dQ m (Q, 0). Hence, the existence of the Chern-Simons form is by itself a proof that C m is 
an exact differential. Moreover, this means that B? can be thought of as providing the local coordinate description 
of the Chern-Simons form. In Appendix [El we show how L m can be locally rewritten as £ m (0) = dQ{Vl, 0): 



c m = (dn AiBi + n\ a si SiBi ) a • • • a (dn A ™ B ™ + n A s ™ a q, s ™ b A a e AlBl ... AmBm 

,1 ' (25) 

= md / dtn aibi A O a t 2b2 A • • • A O a t mbm A e aifcl ... Qmfcm . 



JO 

where t — tdfl t + t 2 fl A SI is the curvature corresponding to the interpolated connection Sl t = tSl which interpolates 
between SI and ST = as t varies between and 1. Although these concepts are somewhat abstract, there exist many 
simple cases where it is actually quite easy to see how these Lagrangians can be rewritten as a total derivative. In 
this sense, it is worthwhile to emphasize certain algebraic aspects of Lanczos-Lovelock Lagrangians, which we will 
now attempt. 

To begin with, because all the indices are contracted with the antisymmetric e tensor, and because the spin 
connections are antisymmetric when working in an orthonormal basis, we have significant freedom in rearranging the 
terms and renaming indices. Hence it is possible to perform a number of calculations without having to explicitly 
worry about the positioning of the indices; Appendix [El is a good example of how advantageous this can be. 

Second, because both dSl and SI A f7 are two-forms, they commute with each other. Combined with the above fact 
regarding the freedom in renaming indices appropriately, this means that we are able to use the binomial theorem to 
rewrite the Lagrangian as, figuratively, 

L m = (dn + si 2 )™ = J2 ( ™ ) {dn) p (n 2 ) m ~ p . (26) 

The p = term in this expansion vanishes identically by antisymmetry. (For example, consider the m = 1 case, where 
the p = term would be 

n A s A fl SB A e AB = 2S7° S A S1 S1 A e i = 0, (27) 

where, in the last equality, we used the antisymmetry of SI and the fact that we are working in two dimensions.) 
Therefore, every term in the binomial expansion will be of the form (dCl) p (£l 2 ) m ~ p with p > 1. 

Third, in 2m dimensions, D'e ai ,,,a 2m vanishes, where D' is the covariant derivative corresponding to any connection 
SI'. Indeed, denoting D g as the covariant derivative built through the metric connection S7 g , we have D g e = 0. This 
means that 

d€ ai ...a 2m = ^g ai e sa 2 ...a 2m + ' ' ' + ^ga 2m £ ai ■ ■ -a 2m -i a 3 

= 0, 

which vanishes because we require all of {s, oi, . . . , aim\ to be different, which is impossible in lm dimensions. This 
then implies D'e = 0. 

Using these ideas, it is straightforward to rewrite some Lanczos-Lovelock Lagrangians as total derivatives directly, 
using the chain rule instead of abstract ideas involving Chern-Simons forms. The m = 2 case is demonstrated in 
Appendix [Dl Since those equations appear, at first glance, challenging, we can here demonstrate the idea on an easier 
example to show that it is actually quite simple: 

dn ab A Sl cd A e abcd 

= d(n ab a n cd ) a e abcd - n ab a dn cd a e abcd 
- d(n ab a n cd ) a e abcd - n cd a dn ab a e abcd 

= ^d(n ab A Sl cd A e abcd ) 
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where in the second-to- last step we equated the second-to- last line with the first, while in the last, we used de a bcd = 0. 

It is therefore clear that one can express C m as an exact differential in the form language. As in the two-dimensional 
case, it is straightforward to use this equation to calculate the set of (non-unique) functions R 3 such that L my J—g = 
djR 3 . One simply has to define an orthonormal basis cj a , use it to calculate the spin connections £l ah , and use those 
in Eq. (l25j) . Therefore, L my f^g can easily be written as a total derivative of a set of functions of the metric and its 
derivatives in any critical dimension. 

A. Spacetimes with a Killing Vector 

If we restrict ourselves to spacetimes which have a Killing vector, then R 3 takes a dramatically simpler form, which 
is worth emphasizing. This fact is a direct consequence of the results obtained in ref.flij and the explicit proof is 
as follows: Suppose the spacetime metric admits a timelike Killing vector £ a which we can choose, without loss of 
generality, to have the components: £ a = (1,0,- •• ,0). Since Q a b = 0, in the critical dimension Q, it follows from 
Eq. ([2} in the case of Lanczos-Lovelock models with X7 a P abcd = that 

L m = 2P 0cde R 0c de (30) 

We now note that 

L m = 2P 0cdb R 0cdb = P 0cdb RacdbC = 2P 0cd6 V c V d & = -^=d c (2^—gP 0cdb V d i b ) (31) 

where we used the fact that P abcd is divergence- free on every index and has the same symmetries as the Riemann 
tensor. It follows that L m ^/—g = d c R c with R c — 2y/— gP acdb V d£,b- Using the relations = Tbod and find that 

our final result can be written in the form: L m ^/—g = djR 3 with 

W =-2^—gP a b ^ b . (32) 

In this case W has the interpretation of being a component of the Noether potential corresponding to diffcomorphism 
invariance (see [13] for a discussion of Noether current in the context of action principles.). This is interesting because 
it gives a direct physical interpretation to the local version of the Chern-Simons form in spacetimes which have 
symmetries. 

Obviously, similar results hold for spacetimes which have a spacelike Killing vector; if the Killing vector is taken to 
be, say, £ J = 5\ in a coordinate system, we can obtain a similar expression starting from the identity Q\ = 0. When 
the spacetime has more than one Killing vector, we obtain, in general, more than one choice of R J which is again an 
example of the non- uniqueness pointed out before. The metric in Eq. (|16[) . for example, has three obvious Killing 
vectors and hence one can write simple expressions for R 3 . They will differ from the expression in Eq. ([17]), as well 
from each other by functions f 3 which satisfy djf 3 — 0. Thus, for a large family of spacetimes, we have a very simple 
way to express the Lanczos-Lovelock Lagrangians as total derivatives. 

V. CONCLUSIONS 

We were led to this investigation and writing up of this paper in a fairly detailed, pedagogical, style because of 
significantly different perceptions of this issue amongst researchers. Those who are trained in the differential form 
language and approach general relativity from that perspective find it a rather trivial result that the action for 
Lanczos-Lovelock models is a surface term in the critical dimension. Given the fact that the action is an integral over 
an exact form (see Eq. ([25]) ). they find it rather strange that there should be any controversy in this matter! On the 
other hand, those approaching the subject from more traditional point of view of tensors, metric, Christoffel symbols, 
etc. do not find it so obvious that L^J—g = djR 3 should hold in the critical dimension. This situation is made more 
confusing because we cannot build a vector density R 3 from the metric and its first derivatives. Further, no explicit 
expressions for R 3 seems to be given in the literature (except a couple of incorrect ones which we have discussed), 
contrary to the often expressed hope by our colleagues that such expressions must exist in the literature! All these 
prompted us to present our results in a rather detailed format, making the paper self-contained. 

In order to clarify several contradictory comments in the literature, we have given a complete description of how the 
2m-dimensional Lanczos-Lovelock Lagrangians, and especially the two-dimensional Einstein-Hilbert Lagrangian, can 
be written locally as total derivatives L m ^/—g = djR 3 with R 3 expressible in terms of the metric and its derivatives. 
This R 3 has been calculated explicitly, as a function of the metric and its derivatives, for the two-dimensional Einstein- 
Hilbert case and we have given a constructive algorithm for the general case. We have also shown that R 3 can be 
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thought of as the local version of the Chern-Simons form written in coordinate language and is related to the Noether 
current in spacetimes which have a Killing vector. 
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Appendix A: Review of the Cartan Formalism 



In order to make sure that the equations are simple to follow, and also to establish notation, we shall provide a quick 
summary of the Cartan formalism, based on the language of differential forms, which is used extensively throughout 
this article. This material is available in a number of other sources, such as 0. 

We have a _D-dimensional manifold M with boundary dM. on which we define a set of coordinates { x 3 } . Since 
we will want to introduce a covariant derivative, we attach, locally on every coordinate patch in this manifold, a 
cotangent bundle T* with basis vectors uj a . These basis covectors are related to the coordinate basis {dx 3 } by a set 
of functions lu°-: uj a — uijdx 3 . Next, we define the spin connection one-form fl a b through the first Cartan structural 
equation: 

duj a + n a b A uj b = 0, (Al) 

for a torsion-free derivative. In particular, for a coordinate basis (that is, ut a = Sfdx 1 ), the spin connection fl a b = 
T^dxi is simply made of the Christoffel symbols. This spin connection lets us define a covariant derivative D, whose 
action on (p, q) tensors is given by 

Dv£:% = dv£;; b a » 

+ A V b c «X" + --- + ^A V^- 10 (A2) 

-« C 6xAC*"S SlVC'v 1C ' 

The curvature two-form is then given by the second Cartan structural equation: 

Q a b = dn a b + n a c a n% = l -R a bcd u c a u> d , (A3) 

where R bcd is the Riemann tensor in the w° basis. The Einstein-Hilbert Lagrangian may then be written 

Q ab A *(uj a A oj b ) = 2 y Rt n uW b e ljni ... nD _ 2 



dx rn A dx n A dx ni A ■ • • A dx n °- 2 ( A4 ) 

R^JTV=gd D x 

Ry/^gd D x, 



where e ai ... an is the Levi-Civita symbol. 

One can choose the basis for the cotangent bundle such that it is orthonormal: it is defined as 



9ij = VabUfaj- (A5) 

This formalism encodes all the spacetime information into the functions uif , called vielbeins, while the Minkowski 
metric r) a b is used to raise and lower indices. In particular, if we demand that Dr\ a b = 0, we will find that £l a b 
must be antisymmetric. This means that there will be \D(D — 1) independent spin connections. Therefore, for the 
two-dimensional case which interested us in Section HU we have only one spin connection: fioi- 
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If we define a set of functions /i fcc a by 



Then it is easy to see that we can solve Eq. (|A1[) by choosing 



Indeed, 



(A6) 



(A7) 



1 



a o a c 

-duj a . 



(A8) 



Let's begin by deriving Eq. 



Appendix B: Details of Section III Al 

in detail, using an orthonormal basis such that Qab is antisymmetric: 
Ry/=gd 2 x = <d ab A *{uj a A u b ) 

= (dn ab + n a c a n cb ) a *(w q a u b ) 

= 2 (dn al + n" c A O cl ) A *(uj A w x ) 
= 2dtt 01 A*(w Awi) 

= 2d(O 01 A*(w Awi)) (Bl) 

+ fi 01 Aft s A*(w s Awi) 
+ ft 01 A fii, A *(w A w s ) 
= 2d (Q 01 A*(w Awi)) 

Then, given an arbitrary two-dimensional metric gji, one can define an orthonormal basis for the cotangent bundle 
by choosing the gauge where Uq — 0. The functions wf which solve gtj = ufLJjTfab are then 



-5oi 



a & 
'911 , 

5oo 



(B2) 



1 / 1 

V^.900 

such that we may write relations between the orthonormal basis u' and the coordinate basis dx l : 

lj° = uj^dx + oj^dx 1 lj 1 = uldx 1 



dx u 



dx 1 



(B3) 



Having found our orthonormal basis, we can move on to calculate the quantities fi bc a defined by Eq. (|A6|) . which then 
let us calculate the spin connection f^oi defined by Eq. (|A7[) . Differentiating lo 1 gives us 



duj° 
doj 1 



J o,i, ,0 



u;" A uj 1 



(B4) 
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This means that we have 



Mioo — — 2- 



o, ,1 



0*1 



(B5) 



Moil — — 2- 



which, we must remember, are antisymmetric in the first two indices. Note also that [i a bc = l^ab^sc- Using now 
Eq. (|A7I) and differentiating it, we get 



dn i = [(/ii 00 Wo),i + (/xon^i - Mioo^i),o] dx° A da; 1 . 
Then, because the Einstein-Hilbert integrand can be written as in Eq. (|B1[) . we can write 

RV~9 = doRi + diRi, 



(B6) 



(B7) 



where 



dO o 1 

■tti — f-mo^i ~ Mon^i 



goi 
3oo ' 



3oo,i + 3n,o 



R\ — MlOO^o 



(B8) 



3oi 

— 3oo,o — *3oi,o + 3oo.i 
3oo 



This solution is clearly not symmetric in t and x. This is due to our gauge choice luq = 0; had we instead chosen the 
gauge lu® = 0, we would have found 



Rl = 



3oi „ 

— 3n,i — ^3oi.i + <7n,o 
3ii 



3oi 

3ii,o + 3oo,i 

3ii 



(B9) 



Because R^J—g = djR? is linear, the most general solution can be written as the linear combination of Ri and R2: 



R4=g~ = XdjR{ + (1 - X)djR 3 2 , 



(BIO) 



for a real parameter A, which we can think of as embodying the extra gauge freedom we have in specifying vielbeins. 
The final form of W is then Eq. ©. 



Appendix C: Conformal Transformation Method 



Knowing how R and \f—g transform, and using R^/—g — djR 1 , we can discover how R J transforms under a 
conformal transformation g^ — > g^ = fl 2 gji- Indeed, we have 



RsT^g 



R + -g ef v e v f n - — g ef y e nv f n 



i? + 2V e (3 e/ -V / n 



R^f~g + 2^V e ( 3 e/ ^V/Q 



n 



R^—g + 2j—gd c ( -g ei n % I + 2T a ab ^-g bi ^n 



(CI) 



12 



Hence, W transform as R? -> W = i?- 7 + ^y/^gg Jl Q, ti . 

n a diagonal gauge with 501 = 0, it is straightforward to write 




This suggests, for any gauge, an ansatz of the form: 

R 3 = -^=(Ai9oo,o + A^oi.o + ^3311,0 

V=5 " (C3) 

+ A^oo.i + A|3oi,i + AePu.i). 

The first solution, Eq. (|B8|) of section [Bl was found by setting Wq = 0, which specified our gauge in the tangent space. 
Similarly, we will impose a condition on A^: that it either be a constant, or be a multiple of — ; note that setting 

cjq = had also assumed goo 7^ 0. This appears somewhat arbitrary, and it is, but we are here trying to guess the 
answer and this just so happens to work. Finally, we demand that R 3 have the correct conformal transformation 
properties, as derived in earlier. This yields 

R° 

I? 1 

If we substitute this into (Ry/—g — djR 3 ), which we want to vanish, we will find some leftover terms such as 
[2go 5oi,oi5ii(^l — B)] , which have no other term against which they can cancel. This forces A = B. If we now 
use A — B into our equation, we find Ry^^g — djR 3 = ■ ■] where [. . .] is some function independent of A. This 

sets A = 1 and B = — 1, such that we retrieve the old solution: Eq. (|B8[) . 

To find the second solution, we perform the same steps, but instead of imposing that A^ either be a constant or a 
multiple of we now impose that \\ either be constant or a multiple of After this time setting A = — 1 and 

B = 1, we get the second solution: Eq. (|B9[) . As in Appendix III A[ we combine the solutions R{ and R J 2 into a general 
solution using Eq. (|B10p . This, once again, gives us Eq. © as desired. 



7= ( 511,0 - ^4—Soo,i + - 1)501,1 
'-9 \ 9oo 

-7= ( 900,1 - B— 500,0 + (B - 1)501,0 
'-9 \ 9oo 



(C4) 



Appendix D: Derivation of equation (1151) 

Eq. (|15p can be derived directly by simply rewriting the Lagrangian as a total derivative using the chain rule. Since 
it is enlightening to see how this can be done without making use of the Chern-Simons form, we will provide the 
details of the calculations. Working in an orthonormal frame, the spin connection Q AB will be antisymmetric and the 
Lagrangian can be written 

C = Q ab A <d cd A e a bcd 

= [dn ab a dn cd + dn ab a a n fd + n a e a n eb a dn cd + n a e a n eb a n c f a n fd ] a e abcd . 

Since d 2 = 0, the first term may be written d(Q AB A dQ CD ). The second term can be expanded as 
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bed 



a-f 



= J2 [ d ( nab A fiC / A Qfd ) + Qab A d ( QC e A ^ ed )] A e a , 



lied 



a-f 



= [ d (n ab a a + n ab a dft c e a n ed - n ab a ft c e a <m ed ] a e Q 



a-f 



= Y [ d (tt ab a n c f a + ft ab a <m c a a ft ad + n ab a dft c b a o bd - n ab a ft c a a <m ad - n ab a ft c fc a dn bd ] a e a&cci 



a-/ 



J2 [d(tt ab A fi^ A n fd ) + n ab A (r] aa dn ca A ft ad + W^ cb A n bd - 77 Qa ft CQ A dft ad - r lbb n cb A dft bd )] A e abcd 



a-f 



= ^ i d ( nab A nC f A ^ /d ) - Vaatt ac A dft ba A n ad - rj bb n cb A dft ab A ft 
a-/ 

+ r/ aa ft ad A ft ca A dft ab + T lbb n db A ft cb A dft ba ] A e abcd 
= [ d ( nab A V c f A ft /d ) + dn ab A (r/ aa ft ac A ft ad - T] bb n cb A ft bd + 7? oa ft od A ft CQ - r; bh ft db A ft cb )] A e a . 



bed 



a-f 



= Y [ d (n ab a n c f a ft /d ) - dft Qb a (n c a a ft Qd + n c b a ft bd + ft c Q a n ad + n c b a n bd )] a e a 



/.;('(/ 



a-/ 



= J2 [ d ( nab A ^/ A Qfd ) ~ 2dnab A A ^1 A £ a. 



bid 



a-f 



-d{n ab a n c f a ft /d ) 



a-/ 



A e a b C( j, 



(D2) 

where we used the antisymmetry of e a bcd to ensure that all of a, b, c, d were different, and the antisymmetry of fl ab to 
ensure that / could not be either c or d. The third term is equal to the second by antisymmetrizing the two pairs 
(a o c) and (do d), while the fourth vanishes exactly by antisymmetry of fl ab , as in Eq. (1271) . Noting that the 
covariant derivative of e must vanish, its exterior derivative will be of the form 

dtabed = ^a e sbcd + ■ ■ ■ + ft/fabes- (D3) 

This is identically zero, as antisymmetry requires that s be different than all of a, b, c and d, which is impossible in 
four dimensions. We can now easily write the Gauss-Bonnet Lagrangian as an exact differential: 

C = Q AB A B CD A e ABCD 

= (dn AB + n A E a n EB ) a (dn CD + n c F a q fd ) 

A tABCD 

= (dn AB a dn CD + 2dn AB a n c F a n FD ) a e ABC D 



(D4) 



= d 

= d 



n AB a dn CD + ^i AB a n c F a n FD j a cabcd 



Appendix E: Proof that £ m = dQ 

There exist many derivations of the Chern-Simons form Q; for an in-depth discussion see [l5l - [T7l ]. We will here 
give a slightly different (and more direct) derivation: we want to show that, locally, C rn = dQ. We introduce the 
interpolated connection f2 f = tfi, which interpolates between the connection and the flat connection Q' = 0, as well 
as its corresponding covariant derivative D t and curvature 0t. Then, we define the Chern-Simons form by 

Q = m [ dm aibl A &t 2b2 A • • • A Q a r bm A e aibl ... ambm . (El) 
Jo 
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The Bianchi identity states that D t &t — 0; a quick calculation shows that this is equivalent to the statement that 
A O2) = dfli A f2 2 — fii A df2 2 , and is therefore trivially satisfied: 



D t ef = d{t 2 n a s a o sfc ) + if7 a s a ef + ^ fc s a 



= t 2 d(« a s a o sfc ) + m a s a (t^ sb + t 2 n s t a o tfc ) 

+ to b s a {tdn sa + t 2 n s t a n ta ) 
= t 2 [d(n a s a n sb ) + n a s a d^ sf> + n b s a e«r s ] 
= 0, 



(E2) 



where in the third equality we removed the terms symmetric in (a -H- b), knowing that we will always be contracting 
with e 0l 6 1 ... am (, m . Then, because the covariant derivative of a scalar is equal to its regular derivative, it is clear that 



dQ 



= m^| dtn ai 
Jo 



AOt 2b2 A-A0^- Ae 0l6l ... 0n 



(E3) 



m / dtD t n a ^ A 9^ b2 A • • • A 9^°"* A e 0l6l ... 



Because of the antisymmetry on e, we have fair amount of freedom in reordering the indices in our equation. In order 
to improve the readability of our equations, we shall therefore drop the indices entirely for this derivation. (One needs 
to remember that Vl 2m = 0.) Then: 

dQ = m [ dtD t (fie?- 1 ) = m [ dt{dfl + 2m 2 )(tdfl + ^fl 2 )™- 1 
Jo Jo 

= m / dtt™- 1 V {dQ + 2m 2 ) [ m ){dn) p {m 2 ) m - 1 - p 
Jo p=0 V p J 

r 1 m_1 / _ 1 \ 
= m J dtt™- 1 y m p J {dn) p+1 t m - 1 - p n 2{m - 1 - p) 



P =i 

m— 1 



o^2(m-p) 



m 



.1 m-l / s 



p=0 
ro-2 



(E4) 



m—1 



+ 2mt m - 1 [ ™ + 1 ) (dn) p+1 t m - p - 1 n 2( - m - p - 1) 



Pi '** 

p=l 



m — 1 
P 



m—1 
p+1 



)2(m-p) 



m — 1 
p-1 



m — 1 
P 



(^) P ^ 2(ro " P) 

p=i 



— r 
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where the second-to-last step involves some simple algebra: 



,n I dtt 2m - p - 1 
'o 

1 



(m- 1)! 



(m - 1)! 



(jp — l)!(m — p)\ p\(m — \ — p)\ 
1 



2m — p (p — l)!(m — p — 1)! \_m — p p 

o T7 tt P + 2 (m-p) 

2m — p p\(ra — pj! 

m ! 

p!(m — p)\ 



(E5) 
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